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1 .  Piecewise  Exponential  Estimator  (PEXE) 

We  wish  to  estimate  an  unknown  life  distribution  F  from  censored  data.  More 
specifically,  let  censored  observation  X.  =  min(Y.,Z.);  6 .  =  1  if  Y .  <  Z . ,  0  if 
Y.  >  1^.  (We  assume  throughout  i=  1,  ...,n.)  6^  indicates  whether  the  observation 

is  a  lifelength  (observation  from  F)  or  a  censored  value  (observation  from  unknown 
censoring  distribution  G) . 

The  standard  method  of  estimating  from  censored  data  is  the  Kaplan-Meier 
estimator  (KME)  obtained  by  maximizing  likelihood.  We  have  shown  that  our  PEXE 
estimator  is  superior  to  the  KME. 

(1)  The  PEXE  provides  a  continuous  estimator  of  F  and  a  step  function 
estimator  of  r,  the  failure  rate  function.  The  KME  provides  a  step  function  esti¬ 
mator  of  F  and  no  estimator  of  r.  Most  applied  people  would  feel  more  comfortable 
with  a  continuous  F  estimator. 

(2)  The  PEXE  provides  a  simple  rule  for  estimating  F(t)  for  t  >  maxCXj , . . . , X  ) ; 
the  KME  does  not  provide  any  rule. 

The  two  estimators  are  asymptotically  alike  in  their  behavior.  For  finite 
sample  size,  they  interlace,  showing  that  they  never  deviate  very  much  from  each 
other. 

In  the  near  future,  we  plan  to  investigate  the  piecewise  geometric  estimator 
(PEGE)  for  estimating  discrete  life  distributions.  The  PEGE  results  should 
correspond  closely  to  those  obtained  for  the  PEXE. 
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2 .  Negative  Association. 


One  of  the  extremely  useful  concepts  developed  for  measuring  "relatedness"  * 

is  that  of  (positive)  association.  Developed  originally  for  use  in  reliability, 
it  soon  found  application  in  many  other  areas,  as  described  in  a  very  large  number 
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The  definition  of  negative  association  differs  from  that  of  positive  associa¬ 


tion  in  two  respects. 

Definition.  X  =  (X^  ,X^  , . . .  .X^)  is  negatively  associated  if  CovEf  ()(*' 1  ^  ), g(X^ 2  ^ )  ]  : 
for  all  nondecreasing  functions  f,  g,  and  mutually  exclusive  subsets  X^\  X^2"*  c  x 

(1)  The  pairs  of  subsets  X^1^  ,  X^  are  mutually  exclusive. 

(2)  The  covariance  is  s  0  rather  than  >  0. 

We  expect  that  the  concept  of  negative  association  will  have  almost  as  many 
applications  as  that  of  positive  association  and  will  elicit  a  good  deal  of  further 
research.  Our  initial  paper  [58]  on  the  subject  gives  a  good  start  in  this 
direction . 

5.  h'ifluence  of  Each  Factor  of  the  Bayesian  Model  on  Other  Factors. 

In  [51]  we  show  that  in  the  Bayesian  model  the  larger  the  sample  observations 
are  stochastically,  then  the  larger  stochastically  is  the  posterior  distribution. 

This  type  of  result  can  be  extended  under  appropriate  conditions  to  other 
pairs  of  factors  in  the  Bayesian  model:  the  choice  of  prior  and  the  resulting 
posterior:  the  sample  and  the  likelihood  function;  etc. 

The  resulting  relationships  should  be  of  considerable  value  in  robustness 
studies,  design  of  experiments,  predicting  the  cost  and  accuracy  of  competing 
experiments,  etc. 

4 .  Percentile  Residual  Life  Functions. 

The  mean  residual  life  function  has  been  studied  theoretically  and  applied 


in  practice.  A  closely  related  concept,  the  percentile  residual  life  function 
(PRLF)  ,  received  very  little  attention  until  the  recent  papers  of  Joe  and  Proschan 
f69],  [70],  [78]  and  [80]. 
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Definition.  The  100a  percentile  residual  life  function  of  a  continuous  life 
distribution  F  satisfies: 

F(t  +  lOOa^)  =  lOOa^  for  -“>  <  t  <  00  and  0  <  a  <  1. 

F(t) 

Joe  and  Proschan  first  study  the  basic  probabilistic  properties  of  the  PRLF. 
Then  they  develop  tests  for  the  properties  of  PRLF's.  They  also  produce  estimators 
of  the  PRLF's.  Finally  they  show  how  to  compare  two  life  distributions  on  the 
basis  of  their  PRLF's. 

We  expect  that  the  approach  based  on  PRLF's  will  prove  popular  and  useful  in 
the  analysis  of  data  from  life  distributions  now  that  some  of  the  basic  technology 
has  been  developed.  We  also  expect  further  theoretical  developments  in  this  area. 

5 .  Unified  Theory  of  Monotonicity  in  Selection  Problems. 

In  the  various  selection  problems  treated  and  solved  in  the  past,  it  has  been 
standard  practice  to  demonstrate  by  individual  and  generally  complicated  methods 
that  the  probability  of  correct  selection  (PCS)  is  a  monotone  function  of  the 
various  factors  affecting  PCS.  By  the  use  of  the  recently  developed  decreasing  in 
transposition  (DT)  inequality  theory,  Berger  and  Proschan  developed  a  unified 
method  for  showing  the  monotonicity  of  PCS;  moreover,  the  proofs  involved  were 
very  much  more  simple,  requiring  only  1  or  2  lines  to  obtain  results  of  much 
greater  power  and  generality  than  the  previous  proofs  in  the  literature.  See  [87]. 

6.  Imperfect  Maintenance. 

In  maintenance  models,  it  is  generally  assumed  that  the  maintenance  actions 
required  are  performed  as  scheduled.  Unfortunately,  in  actual  practice,  maintenance 
performance  is  imperfect  due  to  human  error. 


Probabilistic  models  are  now  being  constructed  specifically  incorporating  the 
effects  of  human  error.  The  solutions  of  these  models  are  more  realistic  and 
approximate  the  actual  outcome  in  practice  more  closely  than  previous  models  which 
naively  assumed  no  human  error. 

See  r 75 71  for  a  detailed  mathematical  analysis  in  the  imperfect  repair  model. 
See  [59]  for  a  description  (without  solution)  of  other  imperfect  maintenance  models 


Exact  Results  for  Survival  Function  Estimators. 


Consider  the  randomly  censored  model  where  wc  observe  Z.  -  min(\.,Y), 

1  i  i 

6.  =  I  { X .  <Y.},  where  X.  is  the  time  to  failure  of  the  i*'1  item,  Y.  is  the  time 

lilt  i 

to  censorship,  and  5..  is  1  if  we  observe  a  failure  time,  0  if  wc  observe  a  censored 
time.  Assuming  X  has  distribution  F,  Y  has  distribution  G,  and  X,  Y  independent, 
the  problem  is  to  estimate  F(t)  =  P(X >  t)  when  F,  G  are  completely  unknown.  Related 
problems  are  to  estimate  the  mean  time  or  median  time  to  failure,  in  the  presence 
of  censoring.  A  central  estimator  in  these  problems  is  the  Kaplan-Meier  estimator 
(KME)  which  can  be  expressed  as 


F 

n 


(t) 


nKn(t)  5fi) 

n  C.  UJ  I{Z,  .  >  t} 

i=l  10 


where  Kn(t)  is  the  empirical  distribution  of  the  Z's,  C^n  =  (n  -  i)  (n  -  i +  1) 

Z(l)  <  •••  <  ^  (ri)  are  t*ie  ordered  Z's  and  6^  is  the  6  corresponding  to  Z^-j. 
Virtually  all  research  in  the  past  has  used  the  asymptotic  distribution  of  nl^2^n(0 
to  study  properties  of  the  estimator.  Chen,  Hollander,  and  Langberg  [48]  derived 
exact  expressions  for  the  moments  of  F  (t)  in  the  proportional  hazards  model 

D 

where  G  =  F  for  B  >  0.  This  enables  us  to  study,  under  proportional 


hazards,  the  bias  of  the  KME  and  to  compare  the  exact  variance  of  the  KME  to  its 
asymptotic  variance.  Current  research  uses  the  methods  of  [48]  to  study  exact 
(i.e.,  for  finite  n)  properties  of  competing  estimators  of  F  and  also  estimators 
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of  functions  of  F,  such  as  p(F)  =  /P(x)dx.  These  results  quantify  the  deleterious 
effects  of  censoring  on  the  estimators.  The  techniques  of  [48]  are  also  being 
used  by  other  authors  to  study  exact  properties  of  procedures  for  censored  data. 

8.  Nonparametric  Bayesian  Estimation. 

Ferguson's  Dirichlet  process  prior  enables  the  incorporation  of  prior  infor¬ 
mation  into  the  estimation  of  an  unknown  distribution  F.  If  one  knows  a  priori 
that  the  distribution  being  estimated  possesses  specific  structure  (e.g.,  symmetric 
about  a  point,  or  exchangeable  in  its  coordinates),  then  Ferguson's  Bayes  estimator 
of  F  is  no  longer  Bayes  in  this  "structured"  model,  and  one  can  turn  to  Dalai's 
G-invariant  Dirichlet  to  develop  a  Bayes  estimator  in  the  structured  model.  Hannum 
and  Hollander  [66]  derived  an  explicit  expression  for  the  Bayes  risk  (using  weighted 
squared  error  loss)  of  Dalai's  Bayes  estimator  of  a  symmetric  distribution  under 
a  G-invariant  Dirichlet  process  prior.  They  assess  the  savings  in  risk  attained 
by  incorporating  known  symmetry  structure  in  the  model  and  provide  information 
about  the  robustness  of  Ferguson's  estimator  against  a  prior  for  which  is  is  not 
Bayes . 

9 .  Tests  for  Nonparamet ric  Classes  of  Life  Distributions. 

Most  tests  developed  for  the  IFR,  IFRA,  NBU,  NBUE  classes  have  been  restricted 
to  the  uncensored  case.  Chen,  Hollander,  and  Langberg  [57]  developed  a  test  for 
the  NBU  class  when  the  data  are  censored.  The  NBU  class  arises  naturally  in 
studies  of  replacement  policies,  multistate  systems,  and  shock  models  and  it  is 
important  to  have  a  means  of  assessing  whether  a  distribution  is  NBU  or  not. 

Chen,  Hollander,  and  Langberg  [61]  also  devised  a  test  for  decreasing  mean  residual 
life  alternatives  (DMRL)  for  the  censored-data  case. 
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We  have  also  introduced  new  classes  of  life  distributions  and  corresponding 
tests  for  both  uncensored  and  censored  data  models.  Guess,  Hollander,  and  Proschan 
T 33 H  introduced  the  "increasing  initially,  then  decreasing  mean  residual  life" 
(IDMRL)  class  of  distributions  and  its  dual,  the  "decreasing  initially,  then 
increasing  mean  residual  life"  (DIMRL)  class  of  distributions,  and  showed  how  these 
classes  can  be  used  for  model  building.  Hollander,  Park,  and  Proschan  [26],  [ 89 J 
introduce  the  "new  better  than  used  at  age  tp"  (NBU  -  tp)  class  and  its  dual  the 
"new  worse  than  used  at  age  tp"  (NWU  -  tp)  class.  Here  tp  is  specified.  These 
classes  and  their  corresponding  tests  can  be  used,  for  example,  to  determine  from 
operational  data  whether  an  airplane  engine  after  tp  hours  of  flight  is  stochasti¬ 
cally  as  good  as  a  new  engine. 
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